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Abstract 

The critical behaviour of 3-dimensional disordered systems with magnetic 
field is investigated by analyzing the spectral fluctuations of the energy spec- 
trum. We show that in the thermodynamic limit we have two different 
regimes, one for the metallic side and one for the insulating side with different 
level statistics. The third statistics which occurs only exactly at the critical 
point is independent of the magnetic field. The critical behaviour which is de- 
termined by the symmetry of the system at the critical point should therefore 
be independent of the magnetic field. 

PACS numbers: 71.30.-hh, 05.45.-hb, 64.60.Cn 
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It is now well known that the introduction of disorder into periodic structures has dra- 
matic effects on the properties of the system. In particular 3 dimensional (3D) systems 
exhibit a metal-insulator transition (MIT) as a function of the disorder. A lot of works |I| 
have already been devoted to the study of this phase transition, theoretically, numerically, 
as well as experimentally. While all the methods seem to agree on the fact that the MIT 
should only be characterized by the fundamental symmetries of the system and that one can 
expect a second order phase transition, its description, e.g. with respect to critical disorder 
and critical exponents, still remains a controversial object of discussions 

Theoretically the standard method applied, namely the 2+e expansion, for the calculation 
of the critical exponents turns out to be inadequate for li = 3 0,^. It seems to be a 
characteristic of several theories for disordered systems using perturbative methods, like the 
nonlinear a-model Q or the self-consistent theory 0, to provide quantitative results in the 
weak localization regime (i.e. for small disorder) but to fail in the region of the critical 
disorder. 

Because of this problem numerical investigations in the frame of the Anderson model of 
localization have played an important role in the description of the MIT. Based on extensive 
computations by means of the transfer matrix method (TMM) it is now usually accepted 
that the critical exponent is ~ 1.4 and is independent of the distribution chosen in the 
Hamiltonian [^|[. This result has been confirmed recently 0,0] using a completely dif- 
ferent method, namely the energy level statistics method (ELSM). Meanwhile other TMM 
studies have been carried out on similar models but with magnetic field, which changes the 
symmetry, driving the system from the orthogonal to the unitary universality class, due to 
the break of the time reversal invariance. Surprisingly, in spite of this change of universality 
class, the same value of the critical exponent has been found with and without magnetic 



field and that independent of the strength of the magnetic field [|Tl|,0. In this letter, us- 
ing the ELSM which will again turn out to be very suitable for such a study, we propose 
to explain this surprising observation with symmetry arguments shedding new light on the 
problem. Our conclusions are based on our numerical results for the spacing distribution 



P{s) and the A3 statistics showing that there is a critical enensemble (CE) charateristic for 
the MIT irrespective of the magnetic field which implies that the critical exponent should 
be independent of the magnetic field. 



In order to investigate the MIT with magnetic field we consider the usual Anderson 
Hamiltonian with an additional phase factor in the off-diagonal elements, 

H = J2^^ri\n){n\ + J2 e*^"''"|n)(m|, (1) 

where the sites n are distributed regularly in 3D space, e.g. on a simple cubic lattice, with 
periodic boundary conditions, and 6n^m = ~dm,n = Only interactions with the nearest 
neighbours are considered. The site energy e„ is described by a stochastic variable. In the 
present investigation we use a box distribution with variance o"^ = W'^/12. W represents 
the disorder and is the critical parameter. 9 = describes the case without magnetic field 
with a Hamiltonian invariant under orthogonal transformation, while 9^0 corresponds 
to the magnetic field case which is invariant under unitary transformation. In spite of the 
simplicity of the model, it contains all the relevant properties necessary to describe the MIT 
with magnetic field. A similar Hamiltonian has been proposed by Pandey et al. [|T^] and 
used to study, by ELSM, the transition from the Gaussian orthogonal ensemble (GOE) to 
the Gaussian unitary ensemble (GUE) in a metallic ring pierced by a magnetic flux ||14|| . 

Based on this Hamiltonian, the MIT in presence of a magnetic field will be studied by 
the ELSM, i.e. via the fluctuations of the energy spectrum. This method has already given 
very interesting results in the case ^ = 0, where the MIT corresponds to a transition from 
the GOE to the Poisson ensemble (PE) |[T5|JT6[| which reflects completely uncorrelated energy 
levels in the localized regime. In the thermodynamic limit one obtains two different regimes: 
GOE for W < Wc and PE for W > Wc, which are separated by a critical ensemble (CE) 
TO| , p^ occuring at the critical disorder Wc- For 9^0 one can expect a transition between 



GUE and PE, but the crucial question is what happens in the vicinity of the critical point. 
Before giving the results we shortly review the ELSM. Starting from Eq.(l) the energy 



spectrum was computed by means of the Lanczos algorithm (which is suited to diagonahze 
such very sparse secular matrices) for systems of size M x M x M with M = 13 and 21, 
disorder W ranging from 3 to 80 and phase 6 = O.Itt. The number of different realizations 
of the random site energies e„ was chosen so that about 2 ■ 10^ eigenvalues were obtained 
for every pair of parameters (M, W) which means between 25 and 90 realizations, for which 
the full spectrum has been computed. For the subsequent investigations only half of the 
spectrum around the band center is considered so that the results are not deteriorated by 
the strongly localized states near the band edges. After unfolding the obtained spectrum 
the fluctuations can be appropriately characterized |jT8| by means of the spacing distribution 
P{s) and the Dyson- Metha statistics A3. P{s) measures the level repulsion, it is normalized, 
so is its first moment because the spectrum is unfolded. A3 which measures the spectral 
rigidity is given by 

A3(L) = (ymin T ^^iV(e) - Ae - Bf de\ . (2) 
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where ( means that we average over different parts of the spectrum. 

Using the random matrix theory (RMT) it is possible to calculate P{s) and A3 for the 
two limiting cases of the spectrum, namely the GUE and the PE. For the metallic side one 



obtains 18 
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Pgue{s) c::^ —sexp{ s^) (3) 

TT^ TT 

^'^^^ = h r^^' ~ + ^')SU(r) dr (4) 
where S|.^£;(r) is the variance of the number of levels in a spectral window of width r and 
is given by 
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The formula (3) for Pgue{s) is exact only in the case of 2 x 2 matrices but remains a good 
approximation for the other cases. For A3 there is no analytical solution and the integral 
(4) has to be calculated numerically. 



For the localized case we have 

Ppe{s) = e-% (6) 

^^{L) = ^. (7) 

In Fig.|I| the results for the Dyson-Metha statistics are reported. We find, as expected, the 
GUE and the PE regimes for small and large disorder respectively as well as the continuous 
transition between them as a function of W. In a previous work it was shown for the 



case ^ = that the As-curves are functions of the system size except at the critical point 
where they are size independent. Moreover the curves were shown to move with increasing 
M towards the GOE for W < Wc and towards the PE for W > Wc- In Fig.|^ we observe 
a similar behaviour for 9 = O.lvr, accordingly this time the curves move towards the GUE 
and the PE. The critical value of the disorder where the curves are size independent turns 
out to be Wc — 16.5. Correspondingly, in the thermodynamic limit we expect two different 
regimes with the GUE and the PE, separated by the CE at the MIT. This has now to be 
compared to the CE obtained in the case without magnetic field. In Fig.^ one sees that 
the CE curves are identical and thus the symmetry should be the same, too. Moreover this 
was checked for 6 = 7i/2 which is the "most Hermitian" case and the same results were 
obtained (cp. Fig.^). Only the critical disorder, which is not a universal value, is slightly 
shifted upward towards Wc — 16.65. Such a behaviour has already been noted in a different 
numerical approach ||12|| . 

The CE curve seems to be proportional to L^/^ for small L and becomes linear when 
increasing L. This would mean that the shape of the curve cannot solely be described neither 
by the results of Kravtzov et al. [|1^] nor by the results of Alt'shuler et al. but rather 



would be in agreement with Ref. for small L and with Ref. [ffB] for large L. Moreover it 



has to be stressed again that the CE curve is completely independent of the magnetic field 



in contrast to Ref. 19 



The same results are derived from the study of P{s) in Fig.|^. For small and large 
disorder we again find the GUE and the PE respectively while at the critical point. 



is independent of M and 9. Another interesting point is the controversy pO| , pO| about the 
asymptotic behaviour of P{s). A careful analysis of the shape of P{s) at the critical point 
gives 

P{s) = Asexpi-Bs"") (8) 

with a = 5/4±0.05 and A,B the normalization constants which supports the results obtained 
m Ref. lH]. Moreover one notes that the value of a is in very good agreement with the 



formula a = 1 + l/ud obtained by Kravtzov et al. relating a to the critical exponent u 
and the dimension d of the system. Finally it has to be stressed that the dependence of the 
results in on f3 {(3 = 1 for the GOE and 2 for the GUE) via some coefficients has not 



been derived from the calculation but assumed from the very beginning leaving open the 
question whether those results, at the critical point, depend on the magnetic field or not. 
Figure ^ demonstrates that this is indeed the case. 

As the critical behaviour is determined by the fundamental symmetries of the system 
one should consider the symmetry at the critical point. From the results in Fig.^ and ^ we 
conclude that the magnetic field has no infiuence on the critical behaviour. This means also 
that the critical exponent u = 1.34 ± 0.10, which was obtained for the ^ = case, is 
transferable to the case with magnetic field. 

Finally we compare these results with experiments. Recent measurements performed on 
the persistent photoconductor Alo.3Gao.7As and on uncompensated Si:P p2| suggest. 



although the situation is still not completely clear, that the MIT is effectively independent 
of the magnetic field. About the value of the critical exponent the situation is not really 



clearer but it can be mentioned that in Ref. |22], considering carefully the range of the critical 



behaviour, a critical exponent of ~ 1.3 has been found which would be in agreement with 
our results. This seems to indicate, at least in some cases, that the MIT is in fact driven 
purely by the disorder. Then it is not necessary to include interactions to characterize the 
critical behaviour although these can be important outside the critical regime. It has to be 



noted that some reserves have been made |E3i about the results obtained in Ref. 221 and it 



will be interesting to check carefully this problem in order to clarify the situation. 

In conclusion we have shown, using the statistical properties of the energy spectrum of 
the Anderson model of localization, that the MIT, which is determined by the symmetry of 
the system at the critical point, is not influenced by the magnetic field. This means that 
both cases, with or without magnetic field belong to the same universality class, described 
by the CE, which opens new perspectives in the comprehension of the MIT. Although this 
universality class is defined only for the critical point it is sufficient to fix the properties 



of the critical behaviour. This is in agreement with previous numerical results |lTl| , p!2[| and 



recent experiments [^,^. We point out that these results are valid in 3D systems and do 



not apply to 2D where we know that the MIT is driven by the magnetic field [24 
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FIGURES 

FIG. 1. Dyson- Metha statistics A3(L) for M = 21. A continuous transition from GUE to 
Poisson statistics occurs as a function of the disorder W (denoted on the right-hand side). The 
dotted line shows the GOE result for comparison. 

FIG. 2. Dyson-Metha statistics A3(L) for system sizes M = 13 and 21 and 
W = 14, 16.5, and 19. 

FIG. 3. The Dyson-Metha statistics As{L) at the critical point for different values of 0. 

FIG. 4. Spacing distribution P(s) for M = 21, W = 3,80 and 9 = O.ln. The histograms are 
the numerical results and the full lines reflect the two expected limiting ensembles, namely the 
GUE for the metallic side and the PE for the insulating side. The points o, * represent P{s) at 
the critical point for M = 21 and M = 13. The symbol o denotes the CE for the case without 
magnetic field = 0). 
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